Node-less atomic wave functions, Pauli repulsion and systematic projector 
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A construction of node-less atomic orbitals and energy-dependent, node-reduced partial waves is 
presented, that contains the full information of the atomic eigenstates and that allows to represent 
the scattering properties in a transparent manner. By inverting the defining Schrodinger equation, 
the Pauli repulsion by the core electrons can be represented as effective potential. This construction 
also provides a description of the Pauli repulsion by an environment. Furthermore, the representation 
leads to a new systematic scheme for a projector augmentation. The relation to Slater orbitals is 
discussed. 

PACS numbers: 71.15.-m, 71.15.Dx, 31.15-p 
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I. INTRODUCTION 

First-principles electronic-structure calculations have 
made a large contribution to the understanding of mate- 
rials. Density-functional theory^ has proven itself as an 
extremely versatile and efficient tool to provide quanti- 
tative information of real materials. While a number of 
open challenges such as strong correlations, embedding 
of environments, various property calculations remain, 
much of the methodology has matured and has been im- 
plemented into efficient program packages. 

The development of electronic-structure methods not 
only has the goal of simulating nature, but also aims 
to develop an efficient language that allows to describe 
and rationalize the properties and the behavior of real 
materials in simple terms. This second aspect in turn 
has been an inspiration for the development of algorithms 
that are efficient on the computer. 

Both, the pseudopotential approach^ and the aug- 
mented wave methods^, which have been linked by the 
Projector Augmented Wave (PAW) method^, stand in 
this tradition. The pseudopotential approach emphasizes 
the free-electron-like features, while the augmented wave 
methods provide, for example, a sound basis for a tight- 
binding description of the electronic structure^, which is 
central to the natural language of chemists. 

A central point in discussing the electronic structure 
is the separation of the chemical effects of the valence 
electrons from those of the core electrons. One can say 
that augmented wave methods and the pseudopotential 
approach came to life by providing their own answer to 
this problem. 

The goal of this paper is to provide a tool to exploit 
this separation of core and valence states from a differ- 
ent point of view. This tool is a representation of the 
atomic wave functions in terms of node-less radial func- 
tions. The node-less wave functions provide a description 
of the Pauli repulsion by the core electrons by an effective 
potential and thus provides a canonical pseudopotential. 



The construction described here, is currently used for 
the pseudization of wave functions in the PAW method, 
which is analogous to the construction of pscudopotcn- 
tials. It is also used for the natural choice of tight-binding 
orbitals, which will be used to incorporate strong corre- 
lations into a DFT environment. 

The concept of node-less wave functions have been in- 
troduced by Topiol, Zunger and Ratner— ii2r— to con- 
struct the first ab-initio pseudopotentials. The re- 
sulting pseudopotentials exhibited a divergent term at 
the atomic site which was a disadvantage for plane 
wave calculations. By relaxing the direct connec- 
tion to the atomic wave function, the divergent term 
could be avoided, which lead to the norm-conserving 
pseudopotentials^ used until today. 

In this paper, we extend the concept of Zunger et al. by 
providing a construction that provides direct contact to 
the underlying physics of the atom. We will demonstrate 
how to exploit this concept for the construction of the 
augmentation of the PAW method and how it can be 
used to describe the Pauli repulsion to the core states by 
an effective semi-local pscudo potential. 



II. DEFINITIONS AND MATHEMATICAL 
PROPERTIES 



A. Notation 



Consider the energy- dependent solution \4>{E)) of a 
spherical Schrodingcr equation 



H-E)\(j}{E)) = 



with 



hm = ! 



(1) 



(2) 
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We suppress the angular-momentum indices £, m, but 
assume that the state \<p(E)) is an angular- momentum 
cigenstatc, which can be written as product of a ra- 
dial function R(r, E) and a spherical harmonics Yi^ m (f), 
namely as 



(r\<l>(E))=R(\r\,E)Ye, m (f f ) . 



(3) 



When we impose the condition that the wave function 
vanishes at the boundary of a sphere with radius r x , we 
obtain a series of bound states 



|VO := \</>(E n )) 



(4) 



with energies E n . 

For the time being, we limit the discussion to the non- 
relativistic Schrodinger equation, i.e. 



(r'\H\r) = 5{f-r') 



2m 



V 2 +v(r) 



(5) 



Extensions to the Dirac equation will be presented 
elsewhere.—. 



B. Construction of node-less bound states 

In a first step, we define a discrete series of node-less 
wave function u n {r) for an atom. 

The series is initiated by the lowest bound state for the 
specified angular momentum 

= ■ (6) 

The series is then constructed recursively by 

{H - E n )\u n ) = -\u n -x) (7) 

with the boundary conditions that the functions vanish 
not only at the box radius r x , but also at the origin, i.e. 



lim |r| 



_1 (r|itri} = for n > 1 . 



(8) 



As shown below, the energies in this series are the same 
as the eigenstates of the Hamiltonian and there is an 
explicit transformation between the real wave functions 
\ip n ) and the node-less wave functions \u n ). The trans- 
formation has the form 



E 

m — 1 



.> n & 



E n ) 



(9) 



The product without a term is defined to be equal to one. 
The identity Eq. [5] can be verified by inserting it into the 
Schrodinger equation Eq. [JJ and the boundary condition 
Eq.[2 

Let us inspect the main properties of the node-less 
wave functions. Figure [1] shows the main characteris- 
tics for the 6s wave function of Au in comparison to the 
normal 6s wave function. 
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FIG. 1: Node-less 6s wave functions shown as full lines for 
Gold (above) and its effective potential (below). The dashed 
functions are the full wave function and the full potential. The 
scale on the vertical axis is the potential value in Hartree. 



The radial part of the wave function starts at the origin 
with a high power of \r\. The depression of the wave 
function near the nucleus can be loosely attributed to 
the Pauli repulsion of the core electrons, that repel the 
valence electrons from the core region. 

The leading power of the node-less wave functions near 
the origin is 



(r\u n ) 



2m f 



(2* + 1)1! 



1 



(2^-l + 2n)!!(2n-2)!! 

1 £+2(n-l) e -A„^ m(r0(1 + o( | r ,|2 )) _ (1Q) 

The leading-order term is derived in appendix El The 
exponential function with the yet unspecified parameter 
X n has been added, because it provides a good description 
of the qualitative behavior of the wave function over a 
fairly wide region near the origin. 

Eq. [10] shows that node-reduced partial waves start 
out at the origin with a increasingly high power, as n is 
increased. We found it surprising that no trace of the 
nodal structure is left near the origin. 



C. Construction of node- reduced wave functions 

Now,we introduce energy-dependent partial waves via 
the differential equation 



(H-E)\q„(E)) = -K-i) 



(11) 



We will call them node-reduced partial waves, because 
their number of nodes is n — 1 less than that of the regular 
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partial wave \<p(E)) at the same energy^ This construc- 
tion allows one to effectively remove the core states from 
the shape of the wave function and from the scattering 
properties. 

The boundary conditions at the origin of the node- 
reduced partial waves are chosen analogously to those 
for \u n ), namely 



lim 

|f|->0 



qn(E,r) 
u n {r) 



1 



(12) 



\qi(E)) is identical to the regular energy- dependent par- 
tial wave \4>{E)) with the normalization from Eq. [2] 

Using E = E n and comparing Eq. [11] with Eq. we 
obtain that the bound states of the node-reduced partial 
waves arc the nodc-lcss wave functions defined before, i.e. 



\q n {E n )) = \u r , 



(13) 



The energy-dependent wave function \<fi(E)) can be re- 
constructed from the node-reduced wave function \q n (E)) 
and the node- less bound states \u n ) by 



\q n (E)) + k> n 

ni—1 j=tn 

n-1 



E, 



(14) 



Thus we obtained a partition of the wave function into 
a node-reduced part and a part related to the core wave 
functions, which is again expressed in terms of node-less 
wave functions. 

As shown in appendix [B] the node-reduced wave func- 
tion can be expressed as difference 

\q n+1 (E)) = (\q n (E)) - \q n {E n ))) (15) 

of two node-reduced wave functions \q n (E)) with one ad- 
ditional node. According to Eq. Qj)l we can construct all 
\q n (E)) from the energy-dependent partial wave \<f>(E)), 
if the latter is normalized such that the leading term at 
the origin remains energy-independent. 

Specifically for the bound state energies, the node- 
reduced wave functions are related, as shown in ap- 
pcndix[C] to the energy derivatives of other node-reduced 
wave functions. 



where 



\q n+j (E n )) = ±\q^(E n )) 



ti\E)) :=di\q n (E)) 



(16) 



(17) 



is the j-th energy derivative of \q n (E)). 

Using Eq.[Tni we find that the node-reduced wave func- 
tions are related to each other by a Taylor expansion 
about one bound state energy 



q n (E)) =V \q n+j (E n )) (E-E n y 



3=0 



(18) 



D. Towards a nodal-theorem for node-reduced 
wave functions 



Eq. [JJ] sheds some light on the question of the node- 
less-ness of the wave function \u n ): 

We will show later, that the assumptions of the nodal 
theorem for the node-reduced wave functions is violated 
in the scattering region. Nevertheless, the assumptions 
are valid in the energy region of bound states. It is our 
hope that a variant of the nodal theorem can be found 
that holds for the entire energy region. 

A nodal theorem for \q n (E)) requires one to make the 
assumption that the node- less functions \u m ) with m < n 
do not reach the outer boundary \r\ = r x for the node 
count. This assumption is reasonable for bound states of 
an atom if the outer boundary is chosen sufficiently far 
away. 

With this assumption, the decomposition of the wave 
function in Eq . [T4l shows that \q n (E)) has bound states 
for all energies Ej with j > n. At the lower bound states, 
that is for j < n the node- less wave function \q n (E)) does 
not have a node, but its contribution is removed by the 
energy-dependent factors in Eq. [HI 

Under the same assumption, the regular nodal 
theorem^ shows that nodes of the node-reduced partial 
waves at the outer bound move inward. Thus the number 
of nodes increases by one at each energy Ej with j > n. 

We can also exclude that nodes leave the interval at 
the center, because the leading order of the Taylor ex- 
pansion about the origin is, according to Eq. 1131 energy 
independent. If a node would migrate across the origin, 
the leading order of \q(E)} would change sign in contra- 
diction to Eq. [TBI 

Two additional assumptions need to be made, namely 
that there is no pair- wise creation or annihilation of nodes 
within the interval, and that there are no accidental 
nodes of the node-reduced wave function \q n (E)) at the 
bound state energies Ej with j < n. 

Thus, for reasonable assumptions, these arguments 
show that the bound wave functions are free of nodes. 
In the scattering region at high energies though, the as- 
sumptions made are violated. Numerical calculations 
show violations of this nodal theorem, namely the pair- 
wise formation of two nodes, in the high-energy region, 
which can be seen in Fig. [5] 

As a consequence of the nodal theorem presented here, 
the number of nodes of the node-reduced wave function 
\q n {E) differs from that of the full wave functions by n—l 
for energies above E n and it vanishes below. 
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E. Properties 

1. Slater orbitals 

In 1929 Slater suggested a simple form for node- less 
atomic orbitals^ of the form 

X (r) = Ar n '- 1 c~^ r , (19) 

where n* is an effective main quantum number, Z is the 
atomic number of the atom, s is a screening constant, do 
is the Bohr radius and A is a normalization constant. 

The guiding idea of Slater is related to ours, namely 
to exploit that the number of nodes is of only secondary 
importance for many of the chemical properties, if only 
the outer tail of the wave function is properly described. 
In this regard, our node-less functions can be considered 
an extension of Slater's orbitals. In contrast to Slater's 
orbitals our node- less functions contain, as a set, the com- 
plete information of the wave functions, at the cost that 
they do not have the simple analytical form of Slater's 
orbitals. 

The qualitative behavior and the similarity to Slater's 
orbitals have been investigated by first by Thiemei& and 
Forsfeil. The series of node-less functions is shown in 
Fig.H 

Interesting to see is that the node-less 3d wave function 
of Fe is much closer in space to the 3s and 3p node-less 
wave functions than to the 4s wave function, to which it 
is energetically closer. It has been a general observation 
that the the node-less wave functions group according to 
their main quantum number into sets with similar range, 
and that these sets are well separated from those with a 
different main quantum number. With increasing atomic 
number, the wave functions with same main quantum 
number become even more similar as seen in Fig. |3l 

The Slater orbitals, shown as dashed lines in Fig. [21 
have been parameterized to reproduce the radial density 
maximum, i.e. the maximum of r u„(r), in value and 
position and the logarithmic derivative at a distance 1.5 
times that of the density maximum. The tails of the wave 
function match fairly well, while there are deviations in 
the inner region. Compared to the Slater orbitals, which 
decay with a power time an exponential, the node-less 
wave functions have a simple exponential tail. 



2. Pauli repulsion 

One motivation for the node-less construction was to 
construct parameter-free pseudopotentials to be used for 
example to embed a quantum mechanical calculation into 
an environment described by a simpler theory. By in- 
verting the Schrodinger equation for the node-less wave 
functions, one can construct an effective potential, that 
includes the effect of the Pauli repulsion. For a node- 
reduced wave function, which obeys Eq. II H we obtain 




r 

FIG. 2: Node-less wave functions of the Fe atom and com- 
parison with Slater orbitals as function of the distance from 
the center. Each graph shows in black the orbitals from one 
main quantum number (full black line) and the corresponding 
Slater orbital (orange dashed line). The orbitals of the neigh- 
boring main quantum numbers are show in (green, full lines). 
The wave functions are scaled to the same maximum value. 
The horizontal axis is scaled for each graph individually to 
make the shape similarity evident. The range of radii on the 
horizontal axis extend from the origin to 0.05, 0.35, 1.3, and 
5 ao, respectively from top to bottom. 

the potential with the Pauli repulsion by 

W^^l^f- (20) 
2mu n {r) 

The potential is semi-local, that is a different potential 
acts on each angular momentum channel. Here, only a 
specific angular momentum channel is considered. 

The effective potential can also be constructed from 
the node-reduced partial waves, \q n (E^)*) 7 but this requires 
an energy below the first bound state to avoid poles in 
the potential resulting from the nodes in the wave func- 
tion. The energy derivative of this potential provides an 
overlap operator, that, however, can also be obtained di- 
rectly from the norm of the node-less wave functions as 
compared to the true nodal wave function. 

In Fig. @J the effective potentials are shown for the 
shells with main quantum number 3 to 5 for iron. The 
potentials become more repulsive with increasing angular 
momentum. 

Interestingly, the potentials form groups with similar 
behavior characterized by the main quantum number. 
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FIG. 3: Node-less wave functions for Gold for s,p,d and f elec- 
trons from top to bottom. The wave functions are normalized 
to have identical maximum value 



This implies that a semilocal potential constructed from 
node-less wave functions with the same main quantum 
number is best approximated by a local potential. Such 
choices will also lead to the most stable and transferable 
pscudopotcntials: In the case of Fe, shown in the exam- 
ple, it is advisable to include the 3s and 3p electrons as 
valence electrons. 

Atomic wave functions without nodes, such as 3d wave 
functions, are special, because they do not experience any 
Pauli repulsion, so that the effective potential is equal to 
the total potential. For atoms having wave functions 
without nodes in the valence shell, such as first-row ele- 
ments, 3d transition-metal elements or lanthanides, the 
intrinsic non-locality as judged from our Pauli repulsion 
potentials is strongest. It is these elements for which the 
construction of pscudopotcntials is most difficult. 

Currently, we arc using these effective potentials to 
construct tight-binding orbitals which will be discussed 
in a forthcoming paper—. 



3. Scattering properties 

The shape of the node-less functions makes them suit- 
able helper functions for the construction of pscudo wave 
functions for the PAW method. 

Pseudo wave functions shall be identical to the true 
wave functions in the bonding region. The s- and p- 
wavc function of transition metal ions have their outmost 




FIG. 4: Effective potential including Pauli repulsion for the 
3s3p3d shell (full) the 4s4p4d shell (dashed) and the 5s5p5d 
shell (dotted) of Fes 2 d 6 . For the shell with main quantum 
number 3 and 4, the effective potentials decrease with in- 
creasing angular momentum. For the 3d channel the effective 
potential equals the true potential. 

nodes fairly close to the covalent radius. As a result a 
pseudo wave function follows the all electron wave func- 
tion into the downswing towards the outermost node. 
The requirement of node-less-ness for the pseudo wave 
function requires a sharp turn outside the node and a 
large region with low amplitude. This sharp turn imposes 
large amplitudes in the pscudopotcntial causing stability 
and transferability problems. If the construction is such 
that the pseudo wave function follows the node-less wave 
function this problem is largely avoided. 

A second difficulty of constructing pscudo wave func- 
tion is to avoid artificial core-like states. They can be 
induced by the proximity to semi-core states that are 
excluded from the valence shell. Because the pole of 
the logarithmic wave function of these states may be no- 
ticeable, the pseudized wave functions try to reproduce 
this pole resulting in artificial undesired states. When 
the node-reduced wave functions are used as basis for 
the construction of the pseudo-potential, this problem is 
avoided. 



4- Towards a systematic projector augmentation 

In the PAW method^, the true Kohn-Sham wave func- 
tions I*) are constructed as 

i*>h*>+E(i^>H0«>Wi*) ( 21 ) 

a / 

from auxiliary wave functions |^) by augmenting them 
with the difference of all-electron and auxiliary partial 
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FIG. 5: Phase shift function of Fe as defined by N(E) + \ - 
iarctan((9 r \n(4>(E, r))) as function of energy. N(E) is the 
number of nodes within the radius at which the derivative 
of the wave function is taken. The phase shift functions have 
been displaced individually by integers for better comparison. 
The dotted lines mark integer values. The full lines are the 
phase shifts of the node-reduced partial wave and the dashed 
lines are those partial waves with the original nodal structure. 
From top to bottom, the s,p,d angular momentum channels 
are shown. The energy region between half-integer and the 
next integer value is the band region for this angular momen- 
tum channel. The steps of the dashed lines between -30 and 
-20 H reflect the 2s, 2p shell and the one shortly below zero is 
the 3s, 3p shell. The jumps in the high energy region reflect 
the annihilation of two nodes within the augmentation region. 



waves, \4>a) arid \<Pa) respectively. The coefficients are 
obtained using the projector functions (p a \, that obey a 
biorthogonality condition 



(Pa|<M = <W 



(22) 



with the auxiliary partial waves. 

Node-less wave functions guide us to a more system- 
atic construction of the augmentation of the projector 
augmented wave method. So far, the mapping between 
all-electron and pseudo partial waves has been done for a 
selected grid of energies^. One difficulty is to ensure that 
the choice of auxiliary partial waves is consistent, that 
is to ensure that the augmentation leads to a converging 
series. 

By inspecting some of the conventional pseudization 
schemes, we observed that the difference in shape be- 
tween the resulting energy- dependent auxiliary partial 
waves \ 4> a ) arid the node- reduced wave functions \q n (E a )) 
depends only weakly on energy. This suggests a new ap- 
proach for choosing the projector augmentation. 



We propose the following ansatz for the auxiliary par- 
tial waves, namely 



\4>(E)) = \q n (E)) + \k) 



(23) 



where \k) is energy-independent function that is entirely 
localized in the augmentation region. A convenient way 
to choose \k) is to construct the pseudo partial wave with 
a conventional method for one single energy for each an- 
gular momentum and then to resolve Eq. [53] for \k). 

Following the recipe of the PAW method^, one con- 
structs raw, energy dependent projector functions via 



W{E)) 



■J 
P 

2m e 



v-E)ME)) . 



(24) 



The auxiliary potential v is a local potential, that is iden- 
tical to the true potential beyond the augmentation re- 
gion and smooth inside. The prime has been attached 
to the raw projector functions to distinguish them from 
the final projector functions obtained by enforcing the 
biorthogonalization Eq. [22] 

Insertion of the Ansatz Eq.[53]into Eq. [53] results, with 
the help of Eq. [TT] in 



\p'(E)) = -|u»-i> + )kn(E)) 



2m, 



+ v-E)\k) 



(25) 



Instead of using a grid of energies we can use the Taylor 
expansion of this expression in the energy about some 
value E,,. 



\p' {j) (E„)) 



-\Un-l) + 



P 



2m. 



+ v-E u )\k) 



-\k)5 hl + {v{r)-v{r))\q^(E u )) ,(26) 



5 j,o 



where \p'^'{E v )) is the j-th energy derivative of \p'(E)) 
at the energy E v . Note, that the last term vanishes 
quickly with increasing j, because the difference of the 
potentials is concentrated in the inner region of the atom, 
while the radius, beyond which |g« ') becomes apprecia- 
ble, rapidly shifts further out with each increasing j. 

The Taylor expansion of the auxiliary partial waves 
defined in Eq. [52] is 

\^\E„))=\qU\E v )) + \k}6 jfi . (27) 

Finally, the biorthogonality condition is enforced ei- 
ther by a Gram-Schmidt like bi-orthogonalization or by 
inversion. 

\<p s ) = \^\E V )) 

\ Vj ) = \<t&{E v )) . (28) 

Tests for this class of projector augmentation will be pre- 
sented in a forthcoming paper. 
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III. CONCLUSIONS 

A set of functions has been presented that provides 
the description of the atomic bound and scattering wave 
functions without the obscuring nodal structure. An ex- 
plicit transformation mediates between the original wave 
functions and what we call "node-less" or "node-reduced" 
wave functions. 

This set of wave functions finds applications in the 
construction of pseudopotcntials or the auxiliary partial 
waves in the PAW method. It provides an effective po- 
tential that describes the Pauli-repulsion of the core elec- 
trons without explicit orthogonalization. The node-less 
wave functions will be used to define tight-binding or- 
bitals for the analysis of chemical binding or the calcula- 
tion of correlated materials. 

Acknowledgement: Financial support by the 
Deutsche Forschungsgemeinschaft through FOR 1346 is 
gratefully acknowledged. 



Appendix A: Power-series expansion at the origin 

The behavior of the node-reduced wave functions at 
the origin is analyzed by the power-series expansion of 
the corresponding inhomogeneous differential equation 
Eq. [TT] The radial part of the node-reduced wave func- 
tion is represented as q(E,r) — J^j o-jr e+3 , that of the 
inhomogeneity as w„_i(r) = — J^j 9j r an d the poten- 
tial is expressed as 



v(r) 



Ze 2 



(Al) 



The radial part of the differential equation Eq. [TT] has 
the form 



r' 2 



2m e 
Ze 2 



m=0 



^1 



It results in the following recursion for the coefficients ctj . 



2m„ 



-h 4 (2t+j + i)j 



Ze 2 
47rer 



ij-l 



171 — 



(A2) 



The regular solution of the homogeneous problem 
starts with do, i.e. aj = for j < 0. The irregular 
solution of the homogeneous problem starts with a—2e-i, 
i.e. a,j = for j < —21 — 1. For the inhomogeneous 
differential equation, we can add a solution of the ho- 
mogenous problem, so that the leading term of the solu- 
tion is determined by the inhomogeneity, given that the 



leading power of the inhomogeneity is higher than r . If 
the lowest-order term of the inhomogeneity is g m , the 
lowest-order term of the wave function is of order m + 2, 
i.e. aj = for j < m + 2. 



2m e 

lj = ~h 2 (2£ + j + l)3 



9j-2 



(A3) 



Thus the behavior at the origin is 

n-l 



(2£+l)!l 



/ 2m e 

\W J (2i-l + 2n)!!(2n-2)!! 
x|rf +2 ^Y e>m (f)(l + 0(\r\)) (A4) 



Thus the leading power increases by two with increas- 
ing main quantum number, which describes the repul- 
sion from the origin. The prefactor of the leading power 
is furthermore energy-independent. 



Appendix B: Recursion for energy dependent 
node-reduced wave functions 



Here we show that \q n (E)) for n > 1 can be represented 
as difference between two functions \q n ^i(E)) at different 
energies by Eq. [15] 

We start with a trial function \f(E)) defined by Eq.[T5l 



\f(E)) := (\q„-x(E)) - | g „_ x (£„_!)>) 



E - E n _i 



(Bl) 



and show that it satisfies all conditions for \q n (E)), 
namely differential equation Eqs. [TT] and boundary con- 
ditions Eq. [H 

Using Eq. ETJ together with Eqs. [TT] and [U for the 

node-reduced functions \q n ^i(E)) , we find 



(H-E)\f(E)) = !«„_!) 



(B2) 



which is identical to the differential equation Eq. [TTJ for 
kn(E)}. 

The boundary conditions at the origin are defined by 
the admixture of the homogeneous solutions, that have 
the leading power — I — 1 and £. Because the leading 
power of \q n -i(E)) in r is, due to Eq. IA41 energy inde- 
pendent, and higher or equal to r , the leading power 
of \f{E)) is determined by the inhomogeneity. Thus the 
boundary behavior is determined directly by the differ- 
ential equation, which is identical for \f(E)) and \q n (E)) 

Since \f(E)) and \q n (E)) obey the same inhomoge- 
neous differential equation and the same boundary con- 
ditions, namely value and derivative at the origin, they 
are identical, which proves Eq. 1151 



Appendix C: Energy derivatives of node-reduced 
wave functions 



Here, we show that the energy-derivatives of the node- 
reduced wave functions at the energy E n of a specific 
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bound state are related by Eq. [16J namely 



\q n+j (E n )) = hq^(E n )) 



(CI) 



Term-by-term comparison with the Taylor expansion of 
\q n +i(E)) on the left-hand side yields 



where the q^{E) is the j-th energy derivative of \q n (E)). 

We start with Eq.[l5jfor \q n+ i(E)) and insert the Tay- 
lor expansion for \q n (E)) about E n 

,^ Eg.M \q n (E))-\q n (E n )) 
\q n +i{E)) = — - — 

= f2^\<lL 3 HEn))(E-E n y^(C2) 



\q%(E n )) = -±-\q^){E n )) (C3) 
3 + 1 



Successive application of Eq. IC3I leads to the desired 
result, Eq . I C 1 1 respectively. Eq. [TBI 



Corresponding author: peter.bloechl@tu-clausthal.de 

2 P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

3 W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

4 S. Topiol, A. Zunger, and M. Ratner, Chem. Phys, Lett 
49, 367 (1977). 

5 D. R. Hamann, M. Schliiter, and C. Chiang, Phys. Rev. 
Lett 43, 1494 (1979). 

6 J. Slater, Phys. Rev. 51, 846 (1937). 

7 O. K. Andersen, Phys. Rev. B 12, 3060 (1975). 

8 P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

9 O. K. Andersen and O. Jepsen, Phys. Rev. Lett 53, 2571 
(1984). 

C. F. Melius and W. A. G. Ill, Phys. Rev. A 10, 1528 
(1974). 

1 A. Zunger and M. A. Ratner, Chem. Phys. 30, 423 (1978). 



A. Zunger and M. L. Cohen, Phys. Rev. B 18, 5449 (1978). 
R. Schade, P.E. Blochl, Th. Pruschke in preparation. 
R. Courant and D. Hilbert, Methoden der Mathematischen 
Physik, vol. 1 (Springer, 1924). 
J. C. Slater, Phys. Rev. 36, 57 (1930). 
M. Thieme, Diploma thesis, Clausthal University of Tech- 
nology, 2003. 

C. Forst, Ph.D. thesis, Vienna University of Technology, 
Austria, and Clausthal University of Technology, Germany 
(2004). 

P. Blochl, in preparation. 

This rule is valid for all practical purposes, but it is violated 
in the scattering region. 



